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^ ! Abstract 

We consider a partial trace transformation which maps a multipartite quantum state 
to collection of local density matrices. We call this collection a mean field state. For the 
t— I \ Hilbert spaces (C 2 )®" and C 2 (g) C 2 (g> C 4 the necessary and sufficient conditions under 

which a mean field state is compatible with at least one multipartite pure state are found. 
. Compatibility of mean field states with more general classes of multipartite quantum 

*^ states is discussed. 

1 Introduction and main results 

Consider a system consisting of n quantum particles and denote TiS 1 ' a Hilbert space of states 
>- \ of the z-th particle. A state of the whole system can be most exhaustively described by a vector 
from the Hilbert space 



h = n {1) ® ■ ■ • ® n {n) 

(a pure state). However, sometimes we need not such a detailed description. For example 
if we are interested only in expectation values of one-particle observables, it is sufficient to 
keep track of just n reduced density matrices p^ = tr i ^ 1 (|\I/) ),..., p^ = tr i ^ n (|^ r )(^ r |) 
describing the states of the individual particles (here the subscript i shows the particles to be 
traced out). This is some kind of mean field description because we completely ignore all types 
of correlations (both quantum and classical) between the different particles. 

An arbitrary collection of n density matrices (p^~\ ■ ■ ■ , p*™**) will be referred to as a mean 
field state. Naturally, some mean field states can not be represented as partial traces of any 
pure state. For example let us take n = 2 and H.^ = Ti^ . From the Schmidt decomposition 
we learn that the density matrices p^ and p^ arising from a pure state must share the same 
set of the eigenvalues. If this is not the case, we can say that a mean field state (p^\ p^) is 
not compatible with a pure multipartite state. Conversely, if the eigenvalues of pW and p^ 
match, we can always find a pure state |*) such that p (1) = tr 2 (|tf)(#|) and p (2) =tr 1 (\^){9\). 

*e-mail: serg@cs.caltech.edu 
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The main purpose of this paper is to find the necessary and sufficient conditions for com- 
patibility between mean field states and multipartite pure states in the case of three or more 
particles where Shmidt decomposition can not be applied. For any Hilbert space S denote D(S) 
a set of all Hermitian, non- negative operators with unital trace (a set of density matrices). Now 
we can formulate our first problem. 

Problem 1. Given a multipartite system Tt = Ti^ (g> • • • ® 1iS n \ find a set of all mean field 
states (p^\ . . . , p^), G D(H^), . . . ,p^ G D(H^) which are compatible with at least one 
multipartite pure state |\?) G TC. 

Let us denote the set to be found M pure . Its formal definition is 

M Jmre = {(p« : p^=tr^mm), |*>e«, (*!*> = !}■ 

The compatibility problem may be generalized to the case when a state of the whole system is 
a mixed one with a specified set of the eigenvalues. Denote eig(p) a string of eigenvalues of an 
operator p (to avoid ambiguity, let us agree that the eigenvalues are arranged in the decreasing 
order). Denote also D(S,\) = {p G D(S) : eig(p) = A} a set of density matrices with a 
specified spectrum A. Our second problem is following. 

Problem 2. Given a multipartite system TC = TC^ ® ■ ■ ■ ® T& 1 ^ and a string of eigenvalues X, 
find a set of all mean field states (p«, . . .,p {n) ), p (1) G D(H {1) ), . . . ,p (n) G D(H {n) ) which are 
compatible with at least one multipartite mixed state p G D(TC, A). 
Let us denote the set to be found M(A). Its formal definition is 

M(A) = {(p« ...,p^) : p«=tr^(p), p G D(TC, A)}. 

The Problem 2 can be reduced to the Problem 1 by introducing one additional particle (ancilla) 
with a space of states TC. One suffices to find the set M pure for the extended system and then 
find the subset of M pure corresponding to the fixed eigenvalues of the ancilla. In that sense the 
two problems are equivalent. 



The solution of the Problem 1 for an arbitrary n is known only for the simplest case Tfo' = C 2 
for all i, see 0. This is a system of n spins | or, equivalently, of n qubits. Let (p^\ . . . ,p^) 
be an arbitrary mean field state of the system. Denote A^ the lowest eigenvalue of i-th qubit's 
density matrix pW. By definition, A^ G [0, |]. Then (p^\ . . . , p^ 1 ') G M pure if and only if for 
each % in the range 1,2, ... ,n one has 

A (i) <^A w . (1) 

For instance, if n = 2, one gets: < X^ 2 \ < X^\ that is A^ = A^ 2 ^ as it should be. For 
n = 3 one gets the triangle inequalities: A (1) < A (2) +A (3) , A (2) < A (3) + A (1) , A (3) < A (1) + A (2) . For 
real spin systems, such as magnetics, we have n ^> 1 and the inequalities (JTJ) impose practically 
no restrictions on the mean field states. However these inequalities may be important in 
quantum information theory, in particular in the study of multipartite analogs of the Schmidt 
decomposition PQI2] and polynomial invariants under local unitary transformations [HJHJin|. m 
Section |21 we describe an alternative proof of obtained independently from the work [6J. 

We have solved the Problem 2 only for a bipartite system and only in the simplest case 
Ti^ = = C 2 . This is a system of two qubits. Let us label the qubits by letters A and B. 
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Consider an arbitrary mean field state (p A ,p B ). Denote X A and X B the lowest eigenvalues of 
the density matrices p A and p B . By definition, X A , X B G [0, |]. Let A = {Ai, A2, A3, A4} be an 
arbitrary string of real non-negative numbers, such that Ai > A 2 > A 3 > A 4 and J2t=i = 1- 
In Section El we will show that (p A ,p B ) G M(A) if and only if 

X A > A4 + A3, 

X B > A 4 + A 3 , , s 

X A + X B > 2A 4 + A 3 + A 2 , {Z) 

\X A — X B \ < min{Ai-A 3 ,A2-A 4 }. 

This result also provides a solution of the Problem 1 for n = 3 with TiP^ = Ti^ = C 2 and 

Although we can not attack the Problem 2 in the general case, we can solve its simplified 
version, namely find the convex envelope [8 a of the set M(A): 

Conv(M(X)) = {(pW,...,p^) : p® = tr,^(p), p G Conv(D(H, A))}, 

where Conv(D(7i, A)) is the convex envelope of the set D(7i,X). By definition, any density 
matrix p G Conv(D(TC, A)) can be represented as p = J2 a p a U a poUl for some probabilities 
Pa > 0, = 1) unitary operators ?7 Q and some density matrix p G D(H,X). The set 

Conv(D(7i, A)) is important for example for description of interaction between the system and 
the classical environment. 

The simplified Problem 2 is considered in Section HJ We show that it can be reduced to the 
analogous classical problem concerning multipartite probability distributions. To describe this 
reduction, let us consider as an example a bipartite system of Alice and Bob, Ti. = 7i A ® 7i B . 
Let X A be the j-th element of the string eig(p A ) and X B be the j-th element of the string 
eig(p ) (one can order the eigenvalues by an arbitrary way; this ambiguity will not affect the 
answer). Let d A and d B be the dimensions of the Hilbert spaces 7i A and Ti B . For any bipartite 
probability distribution p AB (i,j) of integer variables i G [l,^" 4 ] and j G [l,<i B ] denote p AB a 
string of d A ■ d B probabilities p AB (l, 1), . . . ,p AB (l, d B ), . . .p AB (d A , 1), . . . ,p AB (d A , d B ) ordered 
by an arbitrary way (this ambiguity will not affect the answer). In Section 0] we will show that 

d B d A 

Conv(M(X)) = {(p A ,p B ) : Xf = ^p^frj), \f = ^P^frj), p AB ~< A}, (3) 

j=i i=i 

where p AB -< X means that the string p AB is majorized by the string A (the majorization relation 
is briefly discussed in SectionHJ). A set of probability distributions p AB (i,j) satisfying p AB -< A 
is known to be a convex one, and a probability distribution p AB {i,j) is an extremal point if and 
only if p AB = A up to some permutation of the elements. This set of probability distributions 
is thus a classical analogue of Conv(D(TC, A)). The generalization of equality (J3J) to the case of 
three or more parties is straightforward. 

2 Mean field states for arrays of qubits. 

In this section we consider a multipartite system consisting of qubits (spins |) and prove 
the equality Let (p^\ . . . ,p < - n - ) ) be an arbitrary mean field state of n qubit system. By 
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appropriate unitary local operators we can always bring the matrices to the standard form 

P« = A«|1><1| + (1-A«)|0)<0|, ie[l,n), (4) 

where A^ is the lowest eigenvalue of p^\ that is A^ G [0, |]. The proof of (JTJ consists of two 
parts. 

(a) Suppose that (p (1) , . . . , p (n) ) G M pure , i.e. p w = tr J y i (|^)(^|) for some pure state |^) G 
(C 2 )®". Let us prove, for example, the inequality \( n > < Yl^Zi^- K n-qubit system is 
described by the pure state the qubits 1, ... ,n — 1 alone are described by the density 
matrix 

p = tr B (|*><*|). (5) 

Taking into account (j3J) we can write 

n— 1 n— 1 n— 1 

Ea w = (*lX)n 1 |>-]|*> = tr(p(X;n 1 [7])), (6) 

where II ! = |1)(1| is a one-qubit projector and the designation Hi[j] means that the operator 
Hi is applied to the j-th qubit. Note that the operator Y^jZi Hi [7] is diagonal in the standard 
qubit basis and counts a number of l's in a basis vector. Its only zero eigenvalue corresponds 
to the state 1 00 . . . 0). For that reason we will introduce an auxiliary operator 

n-1 

= £n 1 [j] + |oo...o>(oo...o| (7) 

Clearly, all eigenvalues of H are greater or equal to unity, i.e. H > I. As the whole n-qubit 
system is described by the pure state |^), the positive eigenvalues of the density matrices p 
and p^ coincide. In particular the highest eigenvalue of p is equal to the highest eigenvalue of 
which is 1 - A (n) . It means that < (1 - X in) )(ip\ip) for any \ip) G (C 2 )®" -1 . Thus 

we conclude that 

n-1 

= tr(pH) - (00 . . . 0|p|00 . . . 0) > tr(p) - (1 - A (n) ) = A (n) , (8) 

j=i 

which is the desired inequality. 

(b) We should prove that any mean field state satisfying the inequalities JTJ does arise from 
some n-qubit pure state. Obviously one and the same mean field state can arise from many 
pure states, so we have some freedom in choosing them. We will consider only even pure states 
(EPS) which are the states of the following form: 

|^) = E c *\ x )' x = (x 1 ...x n ), XiG{0,l}. (9) 

Here B" lien is a set of all binary strings with the even number of l's. The EPS's span the linear 
subspace of dimension 2 n_1 . 

For any choice of the amplitudes c x in @ the reduced density matrices p- 1 ' are diagonal in 
the {|0), |1)} basis and look as 

P (l) = ( E l^l 2 )l°)(°l + ( E W 2 )|1><1|> <e[i,n]. (10) 
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Thus to guarantee that the matrices pW have the standard form (J1J we should only require 
that 

(0|p»|0> > <l|p«|l>, ie[l,n}- (11) 

Besides, the set of all mean field states which arise from the EPS's is a convex one. Indeed, if 
we take an arbitrary pair of mean field states (po , • • • , Po^) an d (pi , • • • , Pi ) such that 

p« =tr, Vi (|*o><*o|) andpf^tr^d*!)^!), z G [1, n], 

for some EPS's 



l^o) = c ^.°l x ) and l^ 1 ) = c z,il x )> 



then for any t G [0, 1] we can realize the linear combination (1— t)(p^\ . . . , pg™' ) )+t(p^, . . . , p^) 
by the EPS 



J,tu euen 

Consequently a set of all mean filed states which arise from EPS's and satisfy the additional 
requirements (JTTj) is also a convex one. 

The mean field state written in the standard form (@J) can be identified with a real vector 
A G R n defined as A = (A^, . . . , A^). The previous statement can be rephrased by saying that 
the vectors A arising from EPS's constitute some convex subset in M. n which we will denote 
Aeps- By definition, < A**) < |, so that Aeps is a convex subset of the n-dimensional cube / 
with the edge length 1/2. Denote V the set of all 2" vertices of J, i.e. the vectors v = (ei, . . . , e n ), 
e G {0, 1/2}. There are exactly n vertices v G V which violate the inequalities ([Tjl. These are 
the vertices containing exactly one 1/2 component. All other 2 n — n vertices satisfy ((!} and 
contain either zero or greater than one 1/2 components. It will be convenient to introduce a 
special designation to distinguish the vertices of the last type: 

n 1 

V* = {veV : v = (e 1 ...e n ), e 4 <^e i; e % G {0, -}, ie[l,n)}. (12) 

We claim that any vector from V* can be realized by an EPS, i.e. V* C Aeps- Indeed, the 
vectors 

(0...0), (iio...o), (iiio...o),...(i...i)ev 

can be realized by the EPS's as follows: 



(0. 


..0) : 




= |o 


...0), 


(Ho. 


..0) : 




= 2" 


•3(|00> + |11)) ® |0...0), 


2 2 2 


..0) : 




= 2~ 


^(OOO) + |011) + |101) + |110>) ® |0...0) 


(I" 


■\) ■ 




= 2" 


2 ExeB« en 



(13) 



All other vectors from V* can be realized by EPS's applying appropriate qubit permutations to 
the states listed above. We know that Aeps is a convex set, so that any vector from the convex 
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envelope ConviV*) also can be realized by an EPS, i.e. Conv(V*) C Aeps- To conclude the 
proof we will need a simple geometrical lemma which claims that the set ConviV*) coincides 
with those part of the cube I which is specified by the inequalities ([Q: 

ConviV*) = {\ El : A W <^A (j) , ie[l,n}}. (14) 

We encourage the reader to verify this statement for n = 2, 3 by drawing a picture. The 
proof of (|T4*|) for an arbitrary n is given in the Appendix A. Thus for any A G / such that 
\w < ^2-,. \v\ i G [l,n] we have proven that A G Conv(V*) C Aeps, so that A can be 
realized by the EPS. It means that any mean field state satisfying can be realized by some 
n-qubit pure state. 

Note that (|T4^) implies also a converse inclusion A E ps C ConviV*) because we already know 
that inequalities (JTJ are necessary for compatibility of the mean field state with a pure state. 
It means that Aeps = Conv(V*). 

3 Mean field states for two qubits 

In this section we consider a bipartite system composed of two qubits and prove that the 
inequalities (J2J) are necessary for compatibility of a mean field state (p A ,p B ) with at least one 
mixed state from D(C 4 , A). A proof that (J2J) is also a sufficient condition for compatibility is 
placed in Appendix B. 

The combinations of Aj's standing in the righthand sides of the inequalities (J2J) can be easily 
found if we consider two classical bits and joint probability distributions of two bits with the 
specified set of probabilities {Ai, A2, A3, A4} which can be arbitrarily assigned to the events 00, 
01, 10, and 11 (this is a classical analog of the set D{C A , A)). Generally there are 4! such 
distributions. In this case the reduced density matrices p A , p B become the partial one-bit 
probability distributions and A A , X B become the minimal values of p A , p B . In these settings 
the inequalities (J2J) are the best estimates on the quantities X A , X B , X A + X B , and \X A — X B \ as 
can be explicitly verified. The nontrivial result is that this estimates also hold in the quantum 
case and the set M(A) is completely specified by them. 

As we will see below, the strict proof of the inequality \X A — A B | < min{Ai — A3, A 2 — A 4 } is 
rather formidable. However in one particular case it can be proved without any calculations. 
Suppose that the minimum here is zero, i.e. Ai = A3 or A2 = A4. Then the spectrum A 
involves an eigenvalue which is degenerated with multiplicity three or four. In this case we can 
represent any density matrix p G D{C 4 , X) as a linear combination of the identity operator and 
a projector on some pure state. Obviously the term proportional to identity operator just shifts 
the spectrums of p A and p B by the same constant. So it doesn't contribute to the combination 
X A — X B . The term proportional to a projector on a pure state produces the reduced density 
matrices p A and p B with the same spectrum according to the Schmidt theorem. So we conclude 
that X A = X B . 

Suppose that (p A ,p B ) G M(A), i.e. p A = tr B (p), p B = tr A (p) for some p G L>(C 4 ,A). We 
should prove that A A , X B , and A satisfy four inequalities (J2J). Let us represent the matrices p A , 
p B as 

P A = x A u A + {i-x A ){i-u A ), P B = A B n B + (i-A B )(/-n B ), (15) 
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where IT 4 , U B are the projectors on the eigenstates of p A , p B with the lowest eigenvalue, so 
that X A = tr(pH [A]), X B = tr(pH B [B]). As usual the qubit acted on by one-qubit operator is 
indicated in the square brackets. 

Suppose now that O is an arbitrary two-qubit hermitian operator with the eigenvalues 
O x > 2 > 3 > 4 . We claim that 

inf tr(Or ] ) = X 1 4 + X 2 3 + X 3 2 + X 4 1 . (16) 

r7eD(C 4 ,A) 

To verify this statement one can consider the variation 5rj = i[h, rj] where h' — his an arbitrary 
infinitesimal hermitian operator. This variation is obviously consistent with the constraint 
77 G D(C 4 ,X). Suppose that the minimum in ()16|) is achieved at the matrix 77* G D(C 4 ,X). 
Then we have 

5tr(07]*) = itr(0[h, 77*]) = itr(h[rj*, O]) = (17) 

for any h. It means that [77*, O] = so that 77* and O are diagonal in the same basis. Therefore 
there is a permutation a G S4 such that 

tr{0 V *) = AiO ct(4) + A 2 CT(3) + X 3 O a{2) + A 4 CT(1) . (18) 

It is easy to verify that the optimal permutation is always a(i) = i and thus the equality (J 16)) 
is proven. 

Let us substitute O = H A [A] into the equality (fTB^) . In this case we have {Oi, 2 , 3 , O4} = 
{1,1,0,0} and thus 

X A = tr(pU A [A]) > inf tr(r]U A [A]) = A 4 + A 3 . (19) 

r?eD(C 4 ,A) 

By the same way substituting O = Il B [B} and O = U A [A] + U B [B] into (USD we get the 
inequalities X B > A 4 + A 3 and X A + X B > 2A 4 + A 2 + A 3 . 
Unfortunately the last inequality of (j2J), that is 

|A A — A B | < min{Ai — A 3 , X 2 — A 4 }, (20) 

can not be proved using the same idea (the analogous method allows to get the estimate 
\X A — X B \ < Ai — A 4 which is too weak). The proof of (|20|) presented below is a bit formidable. 
Let us define polarization vectors for the qubits A and B: 

a a [r]]=tr(a a [A]r]), b a [r]] = tr(a a [B]r]), (21) 

where a x , a y , a z are the Pauli matrices and rj G D(C A ) is an arbitrary density matrix. Denote 
also a = \/a a a a and b = \Jb a b a the absolute magnitudes of the polarization vectors. We 
can express the eigenvalues A A and X B in terms of polarizations as X A = |(1 — a[p\) and 
X B = |(1 — b[p\). Introduce a functional 

F[r ] ] = b[r ] ]-a[r ] ], r] G _D(C 4 ), (22) 

such that \X A — X B \ = ||F[p]|. The inequality (}2T)j) is equivalent to the following one: 

sup F[rj\ < 2 min{Ai - A 3 , A 2 - A 4 } (23) 

77€D(C 4 ,A) 
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Note that an application of the qubit swap transformation changes the sign of F and doesn't 
changes the eigenvalues of rj so the supremum (J2*3j) can not be negative. Let rf G D(C 4 , A) be 
some density matrix at which this supremum is achieved (may be it is not unique). Consider 
the variation 5rj = i[h,rj*] where — h is an arbitrary infinitesimal hermitian operator. This 
variation is consistent with the constraint rj G D(C A , A). We must require that 5F[r]*] = 0. 

Suppose first that both polarization vectors a a = a Q [?]*] and b a = 6 a [r]*] are non-zero. 
Denote a a = a a /a and b = b a /b the normalized polarization vectors. The equation 5F[r/*] = 
then can be written as 

tr(h [77*, a a a a [A] - b a a a [B]\ ) = 0. (24) 
As h may be an arbitrary hermitian operator we conclude that 

a a [a a [A],r ] *]=b a [a a [B], V *]. (25) 

Let us parameterize the density matrix rf as follows: 



V* = \{I + a a a a [A] +b a a a [B] + A a ^a a [A]^[B}), 



(26) 



where A a ^ is some 3x3 real matrix. Using the commutators [a a , a' 3 ] = 2ze a/37 cr 7 we can trans- 
form the equation into 

a?A & ^ = F'A^ s e lSa . (27) 

As the functional F[rj] is invariant under one-qubit unitary transformations, we can assume 
that both qubits are polarized along z-axis, i.e. a — b — (0, 0, 1). Under this assumption one 
can easily check that the most common solution of the equation ()27|) looks as 



A 



s t 
t -s 
r 



(2f 



for some real numbers s, t, r. Substituting it into (J26|) we can represent the matrix rj* in the 
basis {|00), |01), |10), |11)} of C 4 as follows 



1 

4 



1+a+b+r z 

l+a-6-r 

1-a+b-r 

z 1 -a-b+r 



(29) 



where z = 2(s + it). Note that rf is a block diagonal matrix so its eigenvalues {A^, A?i, A3, A4} 
can be easily found: 



A2 



J(l + a-6-r), A^ = ±(1 + r + yj{a + hf + |z| 2 ), 
l(l-a + b-r), 



l(l + r-^{a + b) 2 + \z\ 2 ). 



(30) 



Recall that rj* G D(C 4 , A) so that the eigenvalues A* coincide up to some permutation with 
the eigenvalues Aj. Although this permutation somehow depends upon a, b, r, z (recall that the 
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eigenvalues Aj are ordered in decreasing way) we need not to consider all possible cases. Instead 
we can notice that 

max min-flA^i) - A ct(3) |, |A ct(2) - A ct(4) |} = min{Ai - A 3 , A 2 - A 4 }, (31) 

cr£S 4 

where the maximum is taken over all permutations of four elements. Note also that the lefthand 
side is invariant under the substitution Aj — > A*. Thus we conclude that 

mm{\Xl - X* 2 \, |A* - AI|} < minjAx - A 3) A 2 - A 4 }. (32) 

Taking A* from the table ()30|) we arrive to the inequality 

- min{|a - b\, y/(a + b) 2 + \z\ 2 } < min{Ai - A 3 , A 2 - A 4 }. (33) 

The minimum in the lefthand side is equal to \a — b\ regardless of the parameter z. Consequently 

b- a < 2min{Ai - A 3 , A 2 - A 4 }, (34) 

which proves the inequality (|23|). 

Suppose now that at the supremum of F[rf\ one of the polarization vectors becomes zero. 
The only non-trivial situation is a a [^*] = 0. In this case the requirement 5F[r]*} = under 
the variation 5rj = i[h,r]*] implies that Sa a [77 *] = for all three components a. By analogy 
with ()25|) this requirement implies that 

[a a [A], V *]=0, a = x,y,z. (35) 

It means that 77* commutes with any operator acting only on the qubit A. This is possible only 
if 77* has the form 77* = (J/2) (g) r) B for some one-qubit density matrix i] B . Thus the eigenvalues 
of 77* are 

A 1 = A 2 = ~(l + 6), A3 = A4=i(l-6). (36) 

Note that 77* G D(C 4 , A) only if Ai = A 2 and A 3 = A 4 so the considered situation is not a generic 
one. Nevertheless we have F[r]*] = b = 2(Ai — A 3 ) = 2(A 2 — A 4 ) so that the estimate (|2^j) is 
valid in this case also. 



4 Connection with multipartite probability distributions 

In this section we consider a bipartite system with local Hilbert spaces Tt A and TC B of arbitrary 
dimensions. We show that a set 

Conv(M(X)) = {(p A ,p B ) : p A = tr B (p), p B = tr A (p), p G Conv(D(H, A))} (37) 

can be described in terms of multipartite probability distributions according to (jHJ). The restric- 
tion that the system is bipartite is not essential and is needed only to simplify the designations. 
Using the arguments presented below, one can easily formulate and prove a multipartite ana- 
logue of (jSJ). We start by collecting all relevant facts concerning the majorization relation. 

Consider arbitrary strings of non- negative real numbers p = (pi, . . . ,Pd) and q = (qi, . . . , qa) 
such that Ylt=iPi = Yli=i 9i — 1- Let us arrange the elements of p and q in decreasing order. 
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Denote p\ the z-th element of p from above and q\ the i-th. element of q from above, such that 
Pi > P2 > ' ' ' — Pd an d (Zi > <?2 — ' ' ' — <?d- The string p is said to be majorized by the string 
q, if for each k in the range 1, . . . , d — 1 we have 

k k 

Erf<E?f- (38) 

i=l i=l 

This relation is denoted as p -< q. We will need three well known facts given below: 

(i) p -< q if and only if p^ = J2 aeSd t a q a (i), i — 1, . . . , d, where Sj is the symmetric group, a G 

is a permutation of d elements, t a > and '^ l(7eSd t a = 1. The existence of such decomposition 
is guaranteed by Birkhoff and von Neumann decomposition of a doubly stochastic matrix into 
a convex combination of permutation matrices. 

(ii) Let S be an arbitrary Hilbert space. Consider any pair of density matrices p, 77 G D(S). 
Then eig(p) -< eig{rf) if and only if p = J2 a t a U a r]Ul, where U a are some unitary operators on 
S, t a >0 and J2 a t a = 1, see jZj. 

(iii) Let S be an arbitrary Hilbert space. For any p G D(S) and for any orthonormal basis in S, 
the string of diagonal elements of p is majorized by the eigenvalues of p, i.e. diag(p) -< eig(p). 

Note that according to the property (ii), the requirement p G Conv(D(7i, A)) is equivalent 
to p G D(H) and eig(p) -< A, so we can rewrite (j3*Tj) as 

CW(M(A)) = {(P A ,P B ) : P A = tr B (p), p B = tr A {p), p G D(H), eig(p) < A}. (39) 

Now we are ready to prove the result 0. The proof consists of two parts. 

(a) Let (p A ,p B ) be an arbitrary mean field state such that 

Xt = J2p AB (i,j), Af = £y*(U), (40) 

3=1 i=l 

for some bipartite probability distribution p AB (i,j) satisfying p AB -< A. Consider the eigenvec- 
tors decomposition of the density matrices p A and p B : 

d A d B 

/ = E A ^>^i> p B = J2 x f^)(j B \, (41) 
1=1 j=i 

where \%a) G H a , \]b) G 7~£ B ■ Define a bipartite density matrix p G D{7i) according to 

d A d B 

p = EEp Ae (^')i^^i ® ^>0'ai- (42) 

i=l j=l 

By definition, p A = tr B (p), p B = tr^(p), and eig(p) = p AB -< A, so that (p A ,p B ) G Conv(M(X)), 

see (jnnj). 

(b) Now suppose that (p A ,p B ) G CWw(M(A)), i.e. that p A = ire(p), p B = trA_(p) for some 
p G D(H) such that eig(p) -< A, see fl3~9"j) . Consider the eigenvectors decomposition (}4"T]) and 
define a bipartite probability distribution p AB (i,j) according to 

p AB (z,j) = (^,Jb|p|^,Jb), i e [l,d A ], 3 G [l,d B \. (43) 
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Clearly, this probability distribution satisfies the requirements ([40)1 . By definition, in the basis 
{\iAi3B)}i,i we have p AB = diag(p), so the property (iii) tells us that p AB -< eig(p). Summa- 
rizing, we have: p AB -< eig(p) -< A. Transitivity of majorization relation implies that p AB -< A. 
This completes the proof of the equality (jHJ) . 

For an arbitrary string of eigenvalues A, let us denote P(X) the set of all bipartite proba- 
bility distributions p AB (i,j) such that p AB = A up to some permutation. The set -P(A) is a 
classical analogue of the set D(7i,X). The property (i) tells us that p AB — < A if and only if 
pAB £ Conv(P(X)). Thus the reduction of the Problem 3 to the analogous classical problem 
corresponds to the substitutions (p A ,p B ) — > (X A ,X B ), p — > p AB (i,j), D(TL,X) — > P(X). The 
partial traces corresponds to taking the average over one variable. 

The inclusion (pW, . . . ,p^) E Conv{M(X)) may be regarded as a necessary condition for 
inclusion (p^\ . . . , p^) E M(A). It provides us some knowledge about the sets M(A) and M pure 
for an arbitrary multipartite system. For example consider a system of three particles A,B, and 
C. Let (p A , p B , p c ) be an arbitrary mean field state. The requirements 

(p A ,p B ) E Conv(M(eig(p c ))), 

(p B ,p c ) E Conv{M{eig{p A ))), (44) 
(p c ,p A ) E Conv{M{eig{p B ))) 

are necessary for inclusion (p A , p B , p c ) E M pure . If A,B, and C are the qubits, one can check 
that these requirements are also sufficient. Meanwhile if 7i A = TC B = C 2 and Ti, c = C 4 , 
one can build an explicit example which shows that they are not sufficient. Indeed, take 
eig(p A ) = {0.6,0.4}, eig(p B ) = {0.5,0.5}, and eig(p c ) = {0.3,0.3,0.3,0.1}. One can verify 
that eig(p A <g> p B ) -< eig(p°), eig(p B ® p c ) -< eig(p A ), and eig(p c <g> p A ) -< eig(p B ), so that (jSJ) 
are satisfied. However (p A , p B , p c ) M pure because there is another necessary requirement 
which is violated by this mean field state, see Section El 

5 Conclusions 

We found the necessary and sufficient conditions under which a collection of density matrices (a 
mean field state) can be produced by the partial trace operations from a multipartite pure state 
for the system of n qubits and for the system C 2 ® C 2 <g>C 4 . These conditions may be described 
by linear inequalities on the eigenvalues of local density matrices. We also found the necessary 
and sufficient conditions under which a mean field state can be produced by the partial trace 
operations from a multipartite mixed state whose eigenvalues are majorized by a given real 
vector. 

6 Acknowledgements 

I would like to thank David DiVincenzo and Alexei Kitaev for advices on majorization tech- 
nique. I am also very grateful to Alexei Kitaev for thoroughly reading the manuscript and for 
many suggestions which helped to improve the paper. This work was supported by NWO-Russia 
collaboration program. 



11 



7 Appendix A 



The purpose of this section is to prove the equality (jHJ). Denote I the subset of the cube I 
which is specified by the inequalities ((TJ), i.e. 

I = {A G R n : A w < A ° } ' ^ AW ^2' * G t 1 ' ( 45 ) 

Clearly J is a convex set. The extremal point of I is a vector A G / for which some m > n 
of the inequalities ()45|) specifying / become the equalities and the corresponding system of 
m equations has a unique solution (thus having a rank n). Let us prove that if A G I is an 
extremal point of /, then A G V*. Because any convex set is a convex envelope of its extremal 
points, it will imply that I C ConviV*). 

To each extremal point A G / we can assign a pair of integers (p, q) where p G [0, n] is the 
number of equalities A w = Ei^i A which are satisfied at A and q G [0, n] is the number of 

equalities A w = 0, A w = ~ which are satisfied at A. By definition p + q > n but we can assume 
that p + q = n because at each extremal point there are only n independent equations. Now 
we are ready to list all extremal points of I and check that they all belong to V*. 

(a) The extremal points of the type (p, n — p), p > 3 

We have the equalities A w = E^A^, A^ = E^A^, A« = E^z A(j) for some i^k^l. 
They are compatible with the inequalities A-,- > 0, j G [1, n] only if A = (0 ... 0). 

(b) The extremal points of the type (2, n — 2). 

The equalities A w = Ej^i A > A = Ej^fc A an< ^ ^he inequalities Aj > 0, j G [l,n] are 
compatible only if AW = A^ and A^ = for j i or j ^ k. Thus we have either A w = A^ = 
or AW = \W = i. The extremal points of both types belongs to V*. 

(c) The extremal points of the type (l,n — 1). 

Suppose that AW = Ej^i A fo r some i [1 3 n ]- If we choose the equality A w = then A^ = 
for all j i. If we choose the equality A w = | then we must choose A( fc ) = | for some 7^ z 
and A^ = for j 7M and j ^ k. If we do not choose neither A w =0 nor A w = | then for each 
j ^ i we must choose either A^ = or A^ = |. This choice is compatible with the equality 
AW = £\ \U) only if AW = or A w = I In all these cases we arrive to the extremal points 
which have been already listed in (a) or (b). 

(d) The extremal points of the type (0,n). 

The extremal point of this type is a vertex of the cube / which satisfies the inequalities (Q). 
According to definition (|T2*|) any point of this type belongs to V*. 

So all extremal points of / belong to V* and thus / C Conv(V*). The inclusion ConviV*) C 
/ is obvious because V* d I and J is a convex set. Thus the equality / = ConviV*) is proven. 

8 Appendix B 

In this section we show that if a mean field state of two qubits (p A ,p B ) satisfies the inequali- 
ties (0) for some eigenvalue string A = (Ax, A 2 , A 3 , A 4 ), then it can be represented as p A = tr B (p), 
P B = t r A(p) for some two qubit density matrix p G D(C 4 ,X), i.e. that (p A ,p B ) G M(A). It is 
convenient to draw a region Q defined by the inequalities (J2J) on the plane (A , X B ) for a fixed 
string A, see FIG. 1. We will describe the desired matrix p G _D(C 4 , A) explicitly for each point 
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Figure 1: The region on the plane (\ A ,\ B ) specified by the inequalities for A = 
(0.7, 0.2, 0.1, 0). For an arbitrary A the coordinates of the vertices O, A, B, C, D are: O = (|, |), 
A = (i, \ -min [Ai - A 3 , A 2 - A 4 ]), B = (min [Ai + A 4 , A 2 + A 3 ], A 3 + A 4 ), C = (A 2 + A 4 , A 3 + A 4 ), 
D = |(2A 4 + A 3 + A 2 , 2A 4 + A 3 + A 2 ). 



inside Q. It is sufficient to consider only the case \ A > X B because the constraint p G -D(C 4 , A) 
is invariant under the permutation of the qubits. Also we will assume that Ai — A 3 > A 2 — A 4 . 
The opposite case may be investigated similarly. We will describe three independent families 
of the density matrices p G D(C 4 , A) corresponding to three triangles OCD, OAB, and OBC 
on FIG. 1. Let us define two qubit states 





= T 


-3 [y/l-a |0,0) + v / l + a |M>] , 




= 2" 


"I [VI + a |0,0) - VI - a |1,1)] , 




= 2" 


-§ [Vi-^|o,i) + vi + ^|i,o)] 5 


*4> 


= 2" 


-1 [Vi + 0|o,i)-Vi-0|i,o>]. 



For any parameters a,/3 G [—1,1] this is an orthonormal basis in C 4 , i.e. (\l/j|\l/A = <5jj. 
Consider the density matrix 

p = A x |*i)(*i| + A 2 |^ 4 )(^ 4 | + A 3 |^ 3 >(^ 3 | + A 4 |^ 2 )(^ 2 | G D(C\\). (47) 

The reduced density matrices p A = tr B (p) and p B = tr^(p) have the lowest eigenvalues 

A A = \ {1 - |a(Ai - A 4 ) - (3(\ 2 - A 3 )|} , \ B = l - {1 - |a(Ai - A 4 ) + (3{\ 2 - A 3 )|} . (48) 

One can easily verify that the region of parameters 0<a<l,0</3<ais mapped to the 
triangle OCD on FIG. 1. Next consider the density matrix 

p = Ai|tfi)(#i| + A 2 |^ 4 )(^ 4 | + A 3 |^ 2 )(^ 2 | + A 4 |^ 3 )(^ 3 | G D(C 4 ,A). (49) 
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The reduced density matrices p A = tr^(p) and p B = tr^(p) have the lowest eigenvalues 

X A = \ {1 - HAi - A 3 ) - /3(A 2 - A 4 )|} , \ B = ~ {1 - |a(Ai - A 3 ) + (3(X 2 - A 4 )|} . (50) 

One can easily verify that the region of parameters < /3 < 1, < a < /3 is mapped to 

the triangle OAB on FIG. 1. Now let us define two qubit states 

|* x > = 2-1 [ v ^^|0,0) + v / TT^|l,l)] , 
|tf 2 ) = 2-1 [VT+^ |0o,O) - v 7 !^ |M b 1)] , 

(51) 

|* 3 ) = 2-2 [vT^ 100,1) + v / TT^ |01,O)] , 

|*4) = 2"5 [v^T+a |0, 1) - v 7 !^^ 1 1,0)] , 

where |0 O ) = cos(|)|0) + sin(|)|l) and |0i) = — sin(|)|0) + cos(|)|l). For any parameters 
a G [—1, 1] and tp G [0, 7r] this is an orthonormal basis in C 4 , i.e. (^1^) = 8ij. Consider the 
density matrix 

p = Ai|*i)(*i| + A 2 |§ 4 )(^4| + A 3 |§ 3 )(^ 3 | + A 4 |^ 2 )(^ 2 | G D(C\X). (52) 
The reduced density matrices p A = trs(p) and p B = trA(p) have the lowest eigenvalues 

X A = \{l - M • ^A? 2 + A| 4 + 2A 12 A 34 cos<^ , A B = i [1 - \a\ ■ {\ x + A 2 - A 3 - A 4 )] , (53) 

where Ai 2 = Ai — A 2 and A 34 = A 3 — A 4 . One can easily verify that the region of parameters 
< a < 1, < <£> < 71" is mapped to the triangle OBC on FIG. 1. Our goal is achieved. 
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